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Abstract 

We show that shunt capacitor stabilizes synchronized oscillations in intrinsic Josephson junction stacks 
biased by DC current. This synchronization mechanism has an effect similar to the previously discussed ra- 
diative coupling between junctions, however, it is not defined by the geometry of the stack. It is particularly 
important in crystals with smaller number of junctions, where radiation coupling is week, and is comparable 
with the effect of strong super-radiation in crystal with many junctions. The shunt also helps to enter the 
phase-locked regime in the beginning of oscillations, after switching on the bias current. Shunt may be used 
to tune radiation power, which drops as shunt capacitance increases. 
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Recently THz radiation was obtained from mesa-type layered crystals with intrinsic Josephson 
junctions (IJJ)^ . The number of junctions in that case was not very large, about 600, and power of 
radiation was enhanced by exciting resonance modes inherent to the crystal, which acts as a cavity. 
Part of energy stored in an excited mode leaked outside the crystal as radiation. The limitation of 
such a design is that the radiation frequency is fixed by the crystal resonances and cannot be 
continuously tuned. A general design for high power tunable source of THz radiation based on IJJ 
in layered superconductors was discussed in Refj2, see Fig. 1. The main idea behind this design 
is to get radiation from crystal boundaries from many synchronized IJJ (up to = 10000) biased 
with the DC current. The current induces DC voltage V between neighboring layers and thus 
produces Josephson oscillations with the frequency coj = 2eV/ h tunable by the DC current. It 
was proposed in Ref. 2 to use a crystal with dimensions Ly > c/cuj ^ Lj. to eliminate effect of 
resonance modes in the a;-direction and use metallic screens at \z\ > L2/2 to eliminate destructive 
interference of electromagnetic waves coming from the surfaces x = ±L^/2. When all junctions 
are synchronized, the radiation power emitted from the crystal edge is proportional to A^^ and may 
reach 1 mW from crystal with Ly = 300 fim, Lx = ^ yum and Lz = 40 /im. 

It was shown in Ref. l2l that radiation by the junctions itself can synchronize oscillations. With- 
out radiation, in-phase oscillations in different junctions may be unstable due to excitation of the 
Fiske modes in the layered crystal. Here we propose additional mechanism of synchronization of 
Josephson oscillations by means of external shunt capacitor, see Fig. [TJ It works in a way similar 
to radiation from the crystal, but now all junctions contribute to the electric field inside the shunt 
capacitor. The effect of shunt stabilization of synchronized oscillations in an array of point-like 
Josephson junctions was discussed previously by Chernikov and Schmidt^. Here we generalize 
their results for extended IJJ and find stabilization condition for such systems. We calculate sta- 
bilization effects of both, radiation and shunt, and compare them quantitatively. We show that 
stabilization effects of shunt with moderate capacitance and of radiation in the super-radiation 
regime (large number of junctions, of order 10"^ — 10^) are comparable, while shunt capacitor is 
much more effective in keeping oscillations in different junctions synchronized at smaller A^. We 
demonstrate also that increase of shunt capacitance results in suppression of radiation and thus 
radiation power may be tuned by shunt. 

To account for the effect of external shunt on oscillations in the IJJ, we use the Lagrangian 
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FIG. 1: Stack of intrinsic Josephson junctions shunted by external capacitance. Light green plates are 
metallic screens, superconducting layers are shown by dark green. 



approach formulated in Ref. |4|. The Lagrangian for the system with shunt shown in Fig. [T]is 



ab 

(curlA)^ ^ ^,2 .2 



^(1 - aVl) Vn - ^(1 - cosy^n) - ^Q^, 



drdz ^ ^_ ' + —CsN'if'. (1) 



Here ^n{'^, i) is the gauge-invariant phase difference between the layers n and n + 1, the coordi- 
nates inside layer are r = x, y, phase difference gradients are Vv9„ = {Vx'^n, ^y^n), the London 
penetration lengths are Ac and \ab for currents between layers and inside layers, respectively, ec 
and €ab are high-frequency dielectric constant for electric fields perpendicular to layers (along the 
2-axis) and along layers, i = \ab/s, where s is the interlayer distance, Cq = c/ (^Q^) and Aj = 7s, 
where 7 = Xc/^ab is the anisotropy ratio, and $0 = nhc/e. 

The first two terms account for the electro-chemical and the Josephson energies of the IJJ. 
The factor (1 — aV^)~\ with second discrete derivative V^A„ = + — 2A„ and 

a = e^^s"^(47r)"^(9/i/(9p, originates from the relation between gauge invariant time derivative of 
the phase difference and difference in the chemical potentials, see Ref. jsj: 

^ ~ ~ Pn+l), (2) 



where Vn, Pn = {Ezn — -E^z,n-i)/(4vrs), and yU„ are the potential, the charge density, and the 
chemical potential in the layer n, while E„ is the electric field in the junction n. Eq. ^ results in 
the relation 

E,„(r,t) = (1 - aVl)-\Bd\c/c)^n{r,t), B, = $o/(27rA,f,A,). (3) 

The next terms in the square brackets account for the kinetic energy of the intralayer currents. 
The intralayer current is jn = (c<l'o/87r^A^5)Q„, where Q„ = — V0n — (27r/$o)A„, where we 
introduced the phase 0„ of the superconducting order parameter and the vector potential A„(r) in 
the layer n. The fourth term in the Lagrangian is the energy of magnetic field inside the crystal. 
The term with Qn is omitted because its contribution is negligible at the low frequencies discussed 
here. 

The electric energy inside the shunt capacitor Cg is accounted for by the last term in Lagrangian, 
written in terms of the average phase difference 

The dissipative function is T^jv^n} = TZci^n} + T^abi^n}, where 



Here and aab are the quasiparticle conductivities perpendicular and along the layers, respec- 
tively. 

The Lagrangian and the dissipative function result in the equations of motion for the phases (pn 
and the vector potential A 

d 5C SC 5TI ^ 

at S(j)n OCpn 6(f)n 

and similar equation for A. We write them in the form 

^2 



dr 



[(fn + (3^] + (1 - aV^) (^^c^^ + sin - V„/iy,„ + V^/i^,„^ = 0, 



(8) 



(V^ - r'Tab)hy,n + TabVu^n = 0, (9) 
(V^ - r^Tab)K,n - TabV.^n = 0, (10) 

Here we use reduced coordinates u = x/\j and v = y/\j, reduced time r = tcOp and frequency 
u = ojj/ujp withcjp = c/ (Acy^) as well as reduced magnetic field h = B/i?c and /3 = NCg/Cj. 



Further, Tab = 1 + J^abd/dr, and Cj = 6cLxLy/ {Atis) is the junction capacitance. We introduce 
reduced dissipative parameters Uab = ^T^f^ab/ {I'^^c^p) and = ^t^cTc/ {^c'^p)- The shunt capac- 
itor effectively enhances the capacitance of the Josephson junctions, but only for synchronized 
oscillations; this enhancement is proportional to the number of junctions. 

The typical parameters of optimally doped BSCCO at low temperatures are = 12, s = 15.6 
A, 7 = 500, \ab = 200 nm, the critical current density jc = $oc/(87r^sA^) = 1700 A/cm^, 
i = 130, = 0.2 and Uc = 0.002, The parameter a ~ 0.1 — 1. It was estimated in Ref. |2|that 
with these parameters crystal optimal for radiation should have sizes ^ 40 fxm, Lx^4: — 6 /im 
and Ly > 300 /im. 

The differential equations (flOl) should be completed by the boundary conditions at x = ±L^/2 
and y = ±Ly/2. Using Eqs. (flOl ) and continuity of and By for the time-independent part of the 
phase difference we obtain at these boundaries 

±{2ns/<l>o)By,,. (11) 

The outside magnetic field is created by the bias DC current and by the induced alternating current. 
For Ly ^ Lx we estimate By(x = ±L^j2,y = ±Ly/2) ^ ±2nl /(cLy), while B^ < By. Here 
I = jLxLy is the total interlayer bias current and j is the bias current density. Hence, we estimate 
time-independent phase difference, 

^{y = Ly/2) - ^{y = -Ly/2) < {2Tifs]L^Ly/{c%) ^ 2Tia,ujjL,Ly/c^. (12) 

Here we used the relation j ^ OcEz in the resistive state when voltage is present. The phase 
difference estimated here is very small for crystals with dimensions smaller than cm. Neglecting 
it we use approximation with y-independent phase difference. 

For alternating part of the phase difference, we find boundary conditions by matching electro- 
magnetic fields inside and outside the crystal.- Outside fields, in half spaces |x| > L^/2, obey 
the Maxwell equations, which fix the ratio between electric and magnetic field. Inside the crystal, 
B„(r, t) given by Eq. (flOl) and the electric by Eq. (|3]). When Ly ^ L^,, c/uj, the predominant 
radiation is along the a;-axis. For weak radiation in the y-direction, we can use the boundary con- 
ditions Vt,v9„ = 0. Thus we can omit dependence of (p„ on the y-coordinate also for alternating 
part of the phase differences. Then we obtain the boundary conditions at x = -^Lx/2: 

± hy^nij^) = ^^^^y^^(pmiuj)[\k^\Jo{k^s\n - m\) + ik^NQ{k^s\n - m|)], (13) 
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where we use the Fourier transforms with respect to time, = co/c, and Jo (a;) and No(x) are the 
Bessel functions. 

The following calculations are similar to those in Ref . Q but accounting for shunt contribution. 
We consider high frequencies to ^ 1 and N ^ £ and neglect finite-size effects along the c-axis. 
The equation for y- and n-uniform phase is 

fLx/2 
-U/2 

Here L^. = L^^/Xj. In the limit a; ^ 1 we look for solution 



/Lx/'i 
du(p{u) + + sin - ^V^v^ = 0. (14) 
-1x12 



ip{u,T) = UT + ■r]{u,T), 1] <t: I. (15) 

The equation for rj is 

duii{u) + iy,r]-fVlr] = -sm{ujT), (16) 

Lx 

which should be solved with the boundary conditions at m = ±L^./2 



5.03 
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The solution is 



r] = (l/2)Im[5 + A cos(k^u)] , k^ = u/i (18) 
B = -[io\l + (30 + i(^Jyc\~\ 

A = iC{[\k^ sm(k^Lj2) + iCu cos(k^Lj2)][u{l + (30 + 
^=[l + 2tioC/(k%)]~\ 

where we approximate sin(/c^I/a;/2) ^ k^Lx/2. The first term in r] is the amplitude of synchro- 
nized {y- and n-independent) Josephson oscillations. It drops as (3 increases. The second term 
describes nonuniform electromagnetic wave inside the crystal. It is generated at the boundaries 
due to radiation field. 

To analyze stability of synchronized Josephson oscillations we consider a small perturbation 
9n{u, t) to the solution uniform in n, ip, 

iPn{u,T) = Lp{u,t) + 9niu,T). (19) 



Equations for 9n{u, r) are obtained by linearization of Eqs. (flOl) with respect to 9n{u, r). The term 
cos[r]{r)]9n{u, t) in the linearized equation couples harmonics with small frequency Q to high- 
frequency terms at f2 ± w. At w ^ 1 we can neglect coupling to the higher frequency harmonics 
± mco with m > 1 and represent the phase perturbation as 



9q + ^ 9q± exp(±iu;r) 



cos(gn)e 



where q = nk/N, k = 1,2, ...N. The complex eigenfrequencies are assumed to be small, 
<^ to. We will find them and also conditions when Im[fi] < for all q (stability condi- 
tion). Substituting 9n into linearized equations (f8l)-(fT0l). excluding oscillating magnetic fields and 
separating the fast and slow parts, we obtain for q ^ the coupled equations 



1 + aq 

1 + aa 



+ ivc^ — C 



9a + Gq V„6'g - 



9q± + G„ +V„9, 



2 T-r2/ 



9a 



g,± V„C7g± 



C = {cosri)r ^ Re[r/^]/2 ^ -(l/2)Re 



a, = 2a(l - cosg), Gg^^ ~ ^ 



2(1 — cos q) 



+ 



(20) 
(21) 
(22) 
(23) 



i{VL =F U})Vah 

Here Gq = Gq^. The boundary conditions for slow and fast components at u = ±Lx/2 and 
q'> -k/N follow from Eq. (fT3l) . 

Finally, we obtain Mathieu equation for slow-varying component with g 7^ 

■ - ' —-^-"^ - " (24) 

(25) 



1 + 



A = Re 



with potential V{u) = Vi{u) + V2{u), 



1+ Ctq 



2[u^ + u^^il + aqY 



V2(u) 



1 



2a;2 
1 

2^ 



Re 
Re 



;i + /30 [/^c. sin(A;^L^/2) + iC^^ cos(A;^L^/2)] 
;i + sin(p+4/2) + K+ cos(p+L,/2)]. 



(26) 



Here = ooGq^ and k+ [{Vl — (X')^G'^^]/[(1 + aq)ecq')]. In the following we consider shunt 
with moderate capacitance /3 < 1. In the lowest order in k^^L^ = uL^ji <^ 1 the part V\{u) 
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reduces to a constant, Vi{u) ^ Re[/C^/ (2w^(l + PC))]^ 



C^{C^ + a) + 1 ecL^ 
+ a)^ + 1 L2 



Treating the coordinate-dependent part of 9g as a small perturbation, we find expression for ^l{q), 

r/3e + /c. W2{q)]\ 



fi^ + iz/cfi ^ <i Re 



1 + /3e ' I + (3^ 



W2{q) = — = = . (29) 

P+L^{p+/k+ + cot(p+L^/2)) 

The first and the second terms in Eq. (l28l) represent the stabilization effects of shunt and radiation, 
and of the long-range interlay er capacitance, respectively. The last term, W2, describes the effect of 
modes Oq± induced inside the crystal due to radiation (parametric excitation of Fiske modes). This 
term leads to instability in the limit of zero dissipation and in the absence of other stabilizing terms. 
It is much smaller than unity. Capacitive shunt and radiation introduce the gap in the spectrum of 
weak distortions and are most effective in stabilization. Their contributions can both reach order 
one for e^L^ < = Ns (in the super-radiation regime) and NCg ~ Cj. As Cj = 60 cm for 
S = 1200 /im^, it is easy to reach this condition. In order to achieve the maximal stabilization 
without sacrificing the radiation power, one needs to choose (3^1. 

We conclude that shunt capacitor stabilizes synchronized oscillations in IJJ stack. The effect is 
particularly useful in crystals with the small junction number or at the initial stages of radiation. 
Shunt may also be used to tune the radiation power. 
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